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r~ | The high kinetic energy outflowing in the jets of microquasars is delivered to the surrounding 

interstellar medium. This energy input can cause the formation of bow shocks and cocoons 
that may be detectable from radio to gamma-ray energies. Evidences or hints of emission from 
jet/medium interactions have been reported for some sources, but little has been done regarding 
the theoretical modelisation of the resulting non-thermal emission. We have developed an 
analytical model based on those successfully applied to extragalactic sources for the interaction 
of AGN jets with their surroundings. Focusing on the adiabatic phase of the growing structures, 
C*~) we give estimations of the expected luminosities through synchrotron, relativistic Bremsstrahlung 

and inverse Compton processes. We conclude that the interaction structures may be detectable at 
C*~) radio wavelengths, while extreme values for the jet kinetic power, the source age and the medium 

fT) density are required to make the emission at high and very high energies detectable. 
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1. Introduction 

About 15 microquasars (MQs) are known so far in the Milky Way [1]. However, this number could 
be much higher if other REXBs do in fact display jets but they are to too faint to have been de- 
tected [2]. MQs show a number of process that resemble those found in extragalactic quasars. For 
instance, in analogy with radio-loud quasar jets impacting on the intracluster medium, one may 
expect strong shocks to develope in the termination regions of MQ jets, from where non-thermal 
radiation should be expected. Hot spots and double-lobe morphologies are common features of the 
powerful extragalactic FR-II sources [3] . In the case of MQs, hints or evidence of such interactions 
have been observed in a number of sources (SS 433 [4], IE 1740.7-2942 [5], XTE J1550-564 [6], 
Cygnus X-3 [7], Cygnus X-l [8], [9], H1743-322 [10], LS I +61 303 [11], Circinus X-l [12]). 
However, only some theoretical work has been done regarding the hydrodynamics of the interac- 
tion, or the resulting non-thermal emission in these systems, as compared with the extragalactic 
case (see, e.g., [13], [14], [15], [16], [17]). 

From a dynamical point of view, MQ jets can transport large amounts of kinetic energy and mo- 
mentum far away from the central binary system. The energy released to the medium can reach 
~ 10 49 — 10 51 erg ([9], [4]). MQ jet termination regions could produce significant fluxes of non- 
thermal radiation even for low radiative efficiencies. 

We have investigated whether a typical MQ fulfills the conditions to be detectable when interacting 
with the surrounding external gas. The main paramenters of the model are the MQ kinetic power 
and age, and the matter densitiy of the environment where the system is embedded in. Furthermore, 
our model requires certain assumptions regarding the magnetic field value in the interaction regions 
as well as the acceleration efficiency at the shock fronts. We have fixed them to reasonable values, 
and we focus our studies on the effects of the dynamical properties of the system and relate them 
to the final radiation output. We have centered our work in the case of a high-mass system in 
order to see the role of a massive and hot companion. In the case of a low-mass system, most of the 
assumptions and conclusions should still be valid, except for the inverse Compton (IC) contribution 
which would be strongly supressed due to the much fainter radiation field of the companion star. 
In next sections we give an overview of the analytical model used to characterize the shocked emit- 
ting zones. The non-thermal fluxes predicted are shown in section Sect. 3. We also provide detailed 
numerical simulations of the hydrodynamical interaction of the jet and the external medium. These 
simulations do ensure the validity of the assumptions taken in the analytical model, although some 
differences arise in the description of the recollimation zone (see Sect. 4). We finally discuss the 
obtained results and give some detectability predictions at different energy bands. 

2. Model description 

The study of the interaction of extragalactic jets with their environments has been extensively 
adressed from the early 70's up to present days. Pioneer works ([18], [19]) stated the main ideas 
regarding the formation and evolution of the luminous radio lobes present in FR-II sources. We 
closely follow these works in adapting these models to galactic MQs interacting with the ISM 
(see also [20]) . We consider a supersonic plasma ejected from the central binary system that is 
decelerated by the external gas accumulated in front of it. This cause the formation of a forward 
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Figure 1: Sketch representing the different interaction zones. The central binary system is located at the 
bottom, and two jets (here, only one jet is shown) emerge from it and extent outwards until they are effec- 
tively decelerated. The jet material that crosses the reverse shock inflates the cocoon, which expands against 
the shocked ISM. A contact discontinuity separates the cocoon and the shocked ISM, the latter being further 
limited by the bow shock. 



shock extending into the ISM, while a reverse shock propagates inwards into the expelled material. 
The shocked plasma travelling in the jet inflates a cocoon after crossing the reverse shock. This 
cocoon is overpressured with respect to the environment although with a much lower mass density, 
and keeps the jet protected against disruption due to external gas entrainment. In this way, most 
of the energy and momentum travelling in the jets are delivered at its termination regions without 
important losses. We concentrate on the adiabatic or Sedov phase [21] of the system, in which the 
jet power is mainly converted in work against the ISM. We take the power and the age of the source, 
as well as the density of the environment, accordingly to keep the interaction structures in this 
phase. For simplicity, we consider a constant environment mass density (Pism) although possible 
medium inhomogeneities should be taken into account in a more accurate model. Furthermore, we 
do not take into account the proper motion of the system. In the case of low-mass systems, the 
proper motions could play an important role regarding the formation of the shock structures when 
this velocity is roughly similar to the bow shock velocity [22]. 

Fig 1. shows the three different interaction zones where non-thermal radiation can be produced. 
The first one is the shell region, which corresponds to the material swept up by the bow shock 
propagating into the ISM. The second one, the cocoon region, corresponds to the shocked material 
of the jet crossing the reverse shock formed in the collision between the supersonic outflow and the 
shell (we note that some mixing with the material present in the shocked shell could be expected 
due to Kelvin-Hemholtz instabilities). The third zone accounts for the recollimation shock formed 
when the lateral pressure of the jet equals the pressure in the cocoon (see also Sect. 4). 
The conical, ballistic jet emerges from the central engine with an opening angle *F ~ r/z ~ 0.1 rad, 
being z the distance to the injection point and r the jet radius. The jet radius becomes approximately 
constant at z con f, the distance where the lateral jet pressure becomes similar to that of the cocoon, 
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P c . The cocoon pressure is balanced with that of the shocked material in the shell, P^,, i.e., P c w 
Pb ~ Gjet x ^mq/H), where Vt, ~ (4/3) %r\ x /b is the bow shock volume, and r\, and Zb are its radius 
and length, respectively. Using the Rankine-Hugoniot conditions for strong shocks, the radius of 
the reconfined jet can be obtained through 



""' x "" T f ltti)M^ ' <2 " 1) 

where y is the adiabatic index of the cocoon material, Vj e t and rj et are the jet velocity and bulk 
Lorentz factor, respectively, and c the speed of light (see, for instance, [20]). After crossing the 
recollimation shock, the jet material still moves at a velocity ~ vj e t- We neglect further recollima- 
tion shocks that may occur (see Sect. 4). 

The length and the velocity of the bow shock at a given source age, fMQ, can be obtained through 
the equations corresponding to the Sedov phase for the growing structures: 

where c\ « 1 is a dimensionless constant that depends on the adiabatic index of the jet material and 
the geometry of the bow shock. The bow shock velocity is: 

Vb = ^ = ^- (2 " 3) 
at j'mq 

The radius of the bow shock, r\,, is found assuming a self-similar relationship between its length 
and width given by R = Zb/Vb, We have taken R = 3 although R can vary from source to source 
[23], [24]. Given the strong compression of the ISM gas, the length and the width of the cocoon 
will be also ~ Zb and ~ Hi, respectively. On the other hand, the size of the reconfinement shock is 
taken as the width of the jet at the location of the recollimation shock, r con f, while its size can be up 
to ~ Zb. We also take r con f to be rougly the size of the accelerator at the reverse shock, while for the 
emitter size we take an extension ~ r^. In the shell zone, the size of the accelerator and the emitter 
are taken both to be ~ r^. Finally, for the sketch of simplicity, we assume homogeneous conditions 
in the cocoon and the shell regions, and we apply a one-zone model to compute the non-thermal 
particle evolution and radiation in these zones. Regarding the magnetic field B in the downstream 
regions, we derive its value taking the magnetic energy density to be ~ 10% (T] = 0.1) of the 
thermal energy density. 

2.1 Non- thermal particle distribution 

A power-law particle distribution N(E) = KE~ P , with an spectral index p ~ 2 (see, for instance, 
[25]), is injected at the reconfinement, reverse and bow shock fronts. The normalization constant 
K is taken such that ~ 1% of the kinetic power flowing in the jets is converted into non-thermal 
energy in the postshock regions right after each shock. On the other hand, the maximum energies 
of the relativistic particles, E max , are calculated equating the energy gain to the different cool- 
ing processes. The particle energy distribution at ?mq is computed taking into account radiative 
(synchrotron, relativistic Bremsstrahlung and IC emission; see, e.g., [26]) and adiabatic losses. Re- 
garding synchrotron losses, the magnetic field energy density considered above is used. Relativistic 
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Table 1: List of the the parameters which remain with a constant value in the analytic model 



Parameter 


Symbol 


Value 


Jet Lorentz factor 


Tjet 


1.25 


Jet half opening angle 




0.1 rad 


Luminosity companion star 




10 39 erg s _1 


Orbital separation 


Rorb 


3 x 10 12 cm 


Self-Similar parameter 


R 


3 


Magnetic equipartition fraction 


n 


0.1 


Non-thermal luminosity fraction 


X 


0.01 



Bremsstrahlung is calculated accounting for the densities in the downstream regions. To compute 
IC losses, we consider only the radiation field from the companion star, an OB star with luminosity 
^star = 10 39 erg s _1 . Adiabatic losses, E = (y/r)E, are computed from the size of the emitters 
and the expansion velocity. At the downstream region after the reconflnement shock there is no 
expansion since the jet radius keeps roughly constant. In the case of the shell region, the expansion 
velocity is vj,, and for the cocoon, it is ~ 3/4 Vb [27]. As noted above, we apply a one-zone model 
to compute the electron energy distribution in the cocoon and the shell. The electron population 
properties, the adiabatic coefficient, and the magnetic and radiation fields are homogeneous in both 
emitting regions. For the electrons injected at the reconflnement shock, their evolution is computed 
assuming that the stellar radiation density decays as <=c 1 /z 1 . Given E max and the evolved electron 
distribution in each emitting region, and accounting for the mentioned radiation mechanisms, we 
can obtain the SEDs for each shocked zone. 

3. Model Results 

The effects on the SEDs when varying the source age, ?mq, the jet kinetic power, Q jeU and the 
ISM density, h\sm are showed in Fig. 2. Emission from the shell (top), cocoon (middle) and jet 
reconflnement (bottom) regions have been computed independently. The contribution of only one 
jet interacting with the ISM is accounted for in the figures. 

Maximum energies have been calculated taking into account the accelerator conditions in each 
zone. In the shell region, E max of electrons range between 2 and 10 TeV, being always limited 
by synchrotron losses. In the cocoon region, itmax can be as high as ~ 100 TeV, due to the high 
velocity of the reverse shock as seen in the jet reference system, and are limited by escape losses, 
while the same mechanism makes the maximum energies in the reconflnement region to be E max ~ 
3-10 TeV. 

The maximum emission levels are obtained through synchrotron emission in the three interaction 
zones, with luminosities that can be as high as ~ 10 32 erg s _1 for powerful sources (right panel in 
Fig. 2). IC emission is the dominant process at the highest energies in the cocoon and reconflnement 
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regions, with luminosities up to ~ 10 29 erg s _1 . Relativistic Bremsstrahlung dominates in the shell 
in this energy range. The spectrum follows the energy distribution of electrons, and can reach 
luminosities ~ 10 31 erg s _1 in the 1 MeV-10 GeV range. This emission channel is on the other 
hand negligible in the cocoon and reconfinement regions, due the low number density of particles 
in these zones. 

The source emission response to the variation of ?mq comes mainly through the larger distances 
from the companion star for older sources and vice-versa. The larger the distance, the fainter the 
companion star radiation energy density u± and the lower the IC contribution. Similarly, higher 
values of «ism make the jet to be braked at lower distances from the central engine. The interaction 
regions are then filled with a higher radiation energy density from the companion star, and the IC 
emission is accordingly enhanced. Finally, all the non-thermal luminosities are proportional to the 
jet power since the total number of radiation particles depend linearly on <2j et . 
Regarding the overall non-thermal emission, the highest radiation output within the set of param- 
eters studied corresponds to the case where <2jet = 10 erg s , fMQ = 10 5 yr and tiism = 1 cm" . 
Taking a source to be at 3 kpc, a radio flux density of ~ 150 mJy at 5 GHz is obtained. Considering 
an angular extension of a few arcminutes, and taking a radio telescope beam size of 10", radio 
emission at a level of ~ 1 mJy beam - 1 can be expected. At the X-ray band, we find a bolometric 
flux in the range 1-10 keV of Fi_iokeV ~2x 10~ 13 erg s _1 cm~ 2 . The X-ray emitting electrons 
have very short timescales and the emitter size cannot be significantly larger than the accelerator 
itself. Although the X-rays from the shell are expected to be quite diluted, the X-rays from the 
cocoon would come from a relatively small region and could be detectable by XMM-Newton and 
Chandra at scales of few arcseconds. In the gamma-ray domain, the flux between 100 MeV and 
100 GeV is Fioo MeV<E<ioo GeV ~ 10~ 14 erg s _1 cm~ 2 , while the integrated flux above 100 GeV is 
^E>iot)GeV ~ 10~ 15 erg s _1 cm~ 2 . These values are too low to be detectable by current Cherenkov 
telescopes. For the weakest jets, i.e, lowest ISM densities and youngest sources adopted in our 
model (<2j et = 10 36 erg s _1 , ?mq = 10 4 yr and wism = 0.1 cm -3 ), the fluxes are strongly suppressed. 
In the radio band, the specific flux is F$ ghz ~ 0.1 mJy beam -1 , and the integrated flux at X-rays 
^i-io keV ~ 10~ 14 erg s" 1 cm -2 , and at gamma-rays FiooMeV<E<iooGeV ~ 2 x 10~ 17 erg s" 1 cm" 2 
and F E >iooGeV ~ 2 x 10~ 16 erg s _1 cm~ 2 . 

4. Hydrodynamical simulations 

In order to check the physical values adopted in the analytical model we have performed numerical 
simulations of the interaction of MQ with the surrounding medium. A two-dimensional finite- 
difference code, which solves the equations of relativistic hydrodynamics written in conservation 
form, has been used (for details, see [29] and [28]). The simulation uses a grid of 320 cells in the 
radial direction and 2400 cells in the axial direction with physical dimensions of 40 x 600 r y An 
expanded grid with 320 cells in the transversal direction, brings the boundary from 40 rj to 500 rj, 
whereas an extended grid in the axial direction, consisting of 440 extra cells, spans the grid axially 
from 600 rj to 900 rj. This enlargement ensures that the boundary conditions are sufficiently far 
from the region of study. The conditions at the boundaries are reflection on the jet axis and in the 
side where the jet is injected, simulating the presence of the counter-jet cocoon, with the exception 
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Figure 2: Obtained SEDs for the shell (top), cocoon (middle) and jet reconfinement (bottom) regions taking 
the values for the jet power Qjet = 10 36 erg s _1 (left panel) and Qjet = 10 37 erg s ~' (right panel), and an 
external gas particle density «ism = 0.1 cm~ 3 (left panel) and «ism = 0.1 cm~ 3 (right panel). Two different 
values for the source age are represented, f M Q = 10 4 yr (solid lines) and 10 5 yr (dashed lines). 
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Figure 3: Lorentz factor (top) and Mach number (bottom) maps resulting from hydrodynamical simulations. 
The simulations were performed using Qj et = 3 x 10 36 erg s _1 , ?mq = 3 x 10 4 yr and «ism=0.3 cm~ 3 . 
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Figure 4: Mass density (top) and pressure (bottom) maps resulting from hydrodynamical simulations. The 
adopted parameters are the same as those of Fig. 3 



of the injection point, where inflow conditions are used. Finally, outflow conditions in the outer 
axial and radial boundaries are used. 

The jet is injected at a distance of 10 18 cm from the compact object, and its initial radius is taken 
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to be ro = 10 17 cm. The time unit of the code is equivalent to 3 x 10 6 s, derived using the 
radius of the jet at injection and the speed of light (tq/c). Both the jet and the ambient medium are 
considered to be formed by a non-relativistic gas with adiabatic exponent T = 5/3. The number 
density in the ambient medium is ni$M = 0.3 cm~ 3 . The velocity of the jet at injection is 0.6c, its 
number density n- s = 1.4 x 10~ 5 cm~ 3 , and temperature T ~ 10 11 K (which corresponds to a sound 
speed ~ 0.1 vjet). These parameter values result in a jet power <2j et = 3 x 10 36 erg s _1 . Figs. 3 and 
4 show the velocity, Mach number, mass density and pressure maps obtained with the numerical 
simulations. 

At the time the simulation is stopped, after evolving « 2.7 x 10 4 yr, the bow-shock is moving at a 
speed 2 — 3 x 10 7 cm s~ l , and has reached a distance ~ 3.6 x 10 19 cm. Initially, the jet expands, 
accelerating and cooling, due to its initial overpressure. When the flow becomes underpressured 
with respect to the ambient medium, a first reconfinement shock is generated close to the injection 
2 x 10 18 cm. The fluid becomes then again overpressured when passing through the shock and this 
process is repeated several times around pressure equilibrium with the external medium. At the 
head of the jet, transonic and subsonic flow velocities result from the increase in temperature and 
decrease in velocity, as the flow crosses the reverse shock. The cocoon and the shell material are 
still in high overpressure with respect to the ambient by the end of the simulation. 
The evolution of the pressure and mass density of the interaction zones in the numerical simulation 
is shown in Fig. 5 as a function of time, for the shell (top) and cocoon (bottom) regions. The 
dotted lines in the figures represent the best fits to these evolution plots, which are used in order 
to extrapolate the values to ?mq = 10 5 yr. The values obtained are in reasonable agreement with 
those found in the analytical model. Pressures range between (2 — 7) x 10~ 10 erg cm~ 3 , while p s h e ii 
stabilizes at ~ 2 x 10~ 25 g cm -3 and p co ~4x 10~ 29 g cm~ 3 for the shell and the cocoon regions, 
respectively. Moreover, we use a self-similar parameter R = 3 in the analytical model. We find 
this value to be in accordance with the results of the numerical simulations, which yield a value 
between 2.5 and 2.7 as can be seen in Fig. 7. Finally, we note that in our model only a strong shock 
at the reconfinement point is assumed, while the hydrodynamical simulations show the existence 
of several conical shocks that develop in the jet when its pressure falls to that of the surrounding 
cocoon (see Fig. 6). Therefore, the non-thermal emission presented in Fig. 2 for the reconfinement 
region should be taken as a rough approximation of the real situation. 

5. Discussion 

Despite we focus on the non-thermal emission from the MQ jet termination regions, thermal 
Bremsstrahlung emission should be expected from the shell. Although the shocks considered here 
are still adiabatic, a non-negligible fraction of the jet kinetic luminosity of up to a few % may be 
radiated via thermal Bremsstrahlung. For bow-shock velocities of few times 10 7 cm s _1 , the ther- 
mal emission would peak at UV-soft X-rays, energy band that is strongly affected by absorption 
in the ISM. In general, observations of the thermal radiation from the interaction structures can be 
still used to extract information of the shell physical conditions (e.g. Cygnus X-l [9]). 
The reason why some MQs show non-thermal emission from the jet/ISM interaction regions, 
whereas in other cases such emission remains undetected, is still unclear. In the context of our 
model, we can study the effects of varying the set of parameters defining the source and their envi- 
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Figure 5: Pressure (left) and mass density (right) evolution in the shell (top) and cocoon (bottom) regions 
as a function of time. A fitting of the results is also shown for both variables (dotted line). This fit serves to 
estimate the simulation values at longer times. This extrapolation is strictly valid only if an homogeneous 
external medium and a constant injection energy rate are assumed. 

ronment, and predict some cases in which the interaction structures may or may not be detectable. 
First of all, the energy input injected to the medium should be high enough, and the jet kinetic 
power varies by several orders of magnitude from source to source. In addition, it could be also the 
case that the density of the surrounding medium is so low that the shell and the cocoon get very 
large and their radiation too diluted to be detectable [30]. Moreover, MQ jets could get disrupted 
at some source age, as it is found in FR I galaxies. If this happened within times much shorter than 
the MQ lifetime, the probability to detect a cocoon/shell structure in the MQ surroundings would 
reduce. On the other hand, some sources may be too far, or the non-thermal fraction too small, to 
detect significant emission from the interaction regions. 

The evolution of the pressure, mass density, the velocities and the Mach number predicted by the 
analytical model are in good agreement with those found in the hydrodynamical simulations for the 
shell and the cocoon regions. Otherwise, several conical shocks may be present within the jet as a 
consequence of pressure adjustments with the surrounding cocoon, instead of the one strong shock 
adopted in the analytical treatment. Finally, the length and width of the structures in the model 
and those found through the numerical simulations are also similar, with a constant ratio R ~ 3 in 
both cases, implying that the physical assumptions used in the analytical treatment are valid to first 
order. 
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Figure 6: Pressure along the jet axis as a function of distance from its injection point, at ~ 10 18 cm, as found 
in the hydrodinamical simulations. Several conical shocks are present, due to the pressure balance with the 
surrounding cocoon: each time the jet pressure falls below that of the cocoon, a shock is formed, keeping 
the jet radius roughly constant until it reaches the reverse shock. 




Figure 7: Evolution of the self-similar parameter R = h/ry, as a function of time. After the pronounced 
initial increase, R remains between 2.5 and 2.7 for most of the simulated time 



The results of this work show that the surroundings of some MQs could be extended non-thermal 
emitters from radio to gamma-rays, although in the VHE regime the fluxes are too low to be de- 
tected. However, taking into account the rough linearity between Q jeU tiism, ?mq ; X an d d~ 2 with 
the gamma-ray fluxes obtained, sources with higher values of these quantities than the ones used 
here (but keeping the adiabatic condition for shocks) may render the MQ jet termination regions 
detectable by present Cherenkov telescopes. In addition, from a comparison with observations, 
the magnetic field and the particle acceleration efficiency in the jet/ISM interaction regions can be 
constrained, giving an insight on the physics of these interaction structures. 

To conclude, it is interesting to note that, although the adopted model is rather simple, it already 
accounts for cases when the sources should remain undetectable and cases in which radiation could 
be detected. 
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